Enhanced vaccine control of epidemics in adaptive networks 
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We study vaccine control for disease spread on an adaptive network modeling disease avoidance 
behavior. Control is implemented by adding Poisson distributed vaccination of susceptibles. We 
show that vaccine control is much more effective in adaptive networks than in static networks due to 
an interaction between the adaptive network rewiring and the vaccine application. Disease extinction 
rates using vaccination are computed, and orders of magnitude less vaccine application is needed to 
drive the disease to extinction in an adaptive network than in a static one. 



Modeling the spread of epidemics on static networks 
is a well developed field, but recent studies have begun 
to account for potentially time-varying network topolo- 
gies. In particular, people may adjust their social be- 
havior in response to the threat of an epidemic. Both 
an SIS (susceptible-infected-susceptible) model [l[ and 
SIRS (susceptible-infected-recovered-susceptible) model 
Q have been studied on an adaptive network in which 
non-infected nodes rewire their links adaptively away 
from infected neighbors and toward other non-infected 
nodes. Such adaptation typically increases the epidemic 
threshold and reduces the number of infectious cases, 
and new disease dynamics and bifurcations are observed. 
Similar results are seen when the nodes reconnect to ran- 
domly selected nodes elsewhere in the network 0, 
Because vaccines are available for many diseases, it is 
desirable to examine the interplay between adaptively 
fluctuating social contacts and vaccination of susceptible 
individuals. 

Almost all diseases exhibit randomness resulting in ob- 
served fluctuations, as in 

[1 S, 0, !, i, El- As diseases 
evolve in large populations, there is the possibility of fi- 
nite time extinction 0, El, d El Q Gil • Extinction 
occurs where the number of infectives become so small 
that there is insufficient transmission to keep the disease 
in its endemic state [H, U3, EH ■ Fluctuations cause the 
disease free equilibrium (DFE) to be reached in a finite 
time. Such an extinction process occurs even when the 
DFE is unstable. Populations based upon adaptive net- 
works further complicate the problem, since social dy- 
namical situations, such as disease avoidance strategies, 
can cause the endemic and DFE to be bistable 0]. 

A major characteristic of fluctuation-induced extinc- 
tion in stochastic models for globally connected large 
populations is the extinction rate. Viewing disease fade- 
out as coming from systems far from thermal equilibrium, 
finite population extinction rate laws have been derived 
m SIS [H[il and SIR [2(| models. In contrast to vaccine 
strategies which stabilize the DFE [H [H, [H, [H [H, HH , 
periodic pulsed vaccination was generalized to a random 
Poisson strategy, which exponentially enhances the rate 
of extinction 1271. 



Other vaccine strategies have been examined in a vari- 
ety of static network geometries. Targeting of high degree 
nodes is widely recognized as more effective than ran- 
dom vaccination, including in scale free networks (28|.|29|. 
small world networks [301 ] . and more realistic social net- 
work geometries [3l|. Because targeting the highest de- 
gree nodes requires full knowledge of the network geome- 
try, other strategies based on local knowledge have been 
developed. For example, vaccinating a random acquain- 
tance of a randomly selected node tends to favor high 
degree nodes and is again more effective than random 
vaccination [33 |. When vaccine is very limited, outbreaks 
can be minimized by fragmenting the network via a graph 
partitioning strategy which requires less vaccine than tar- 
geting high degree nodes [331 ] . Here we will use a random 
vaccination strategy and find that in conjunction with 
adaptive rewiring, it is extremely effective. 

We create a new model by modifying an SIS model [l[ , 
adding a vaccine class (V). Individuals are connected in 
a network, with non-infected nodes rewiring adaptively 
to reduce connections with infected nodes. 

The transition probabilities are as follows. A suscepti- 
ble node becomes infected with ratepiV/ jn b r , where Ni^,- 
is its number of infected neighbors. An infected node re- 
covers (to susceptible) with rate r. Vaccination occurs 
in Poisson-distributed pulses with frequency v, and each 
pulse vaccinates a fraction A of the susceptible nodes 
[341 ]. We assume that the vaccine is not permanent, so a 
vaccinated node becomes susceptible again with rate q, 
the resusceptibility rate. 

While the epidemic spreads, the network is rewired 
adaptively. Because the vaccine wears off, we assume 
that vaccinated nodes are uncertain of their infection sta- 
tus and thus rewire in the same way as susceptibles. If a 
link connects a non-infected node to an infected node, 
that link is rewired with rate w to connect the non- 
infected node to another randomly selected non-infected 
node. Self links and multiple links between nodes are 
excluded. 

We simulate the full adaptive network via Monte Carlo 
simulation in a similar fashion to Q • Vaccination events 
occur with average rate v, and in each event AN$ suscep- 



2 



tibles are selected randomly for vaccination, where Ns is 
the total number of susceptibles. Results are presented 
here for N = 10 4 nodes and K = 10 5 links. 

We have developed a mean field theory for the dynam- 
ics of nodes and links following [l|, Q • Pa denotes the 
probability of a node to be in state A, where A is either 
S, I, or V. Pab denotes the probability that a randomly 
selected link connects a node in state A to a node in 
state B. T] = T](t) is the Poisson distributed vaccination 
with mean frequency v and amplitude A. We obtain the 
following set of stochastic differential equations for the 
nodes: 
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The mean field equations are analyzed as a stochastic 
system with random r)(t) as above, as well as for a deter- 
ministic system using the mean vaccine rate (f)(t)) = vA. 

Sample time series are shown in Figure [1] for a run 
in which the epidemic became extinct. Because we use 
proportional vaccination, as vaccination begins to lower 
the number of infectives and the number of susceptibles 
increases, subsequent vaccine pulses vaccinate a larger 
number of nodes. 

We studied the longtime behavior for the case of longer 
lived endemic states. The deterministic mean field model 
predicts stable steady state dynamics for the static net- 
work, even when vaccine is applied. Small regions of 
oscillatory behavior have been predicted for the adap- 
tive network without vaccination [![, and these regions 
are predicted to be much larger when vaccination is ap- 
plied. For the resusceptibility q used here, oscillations 
are generic 
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2(d) compare the dependence of the mean 



infective level on the vaccine frequency for static and 
adaptive networks [36|. For the full static (Erdos-Renyi) 



network system (Fig. 2(b) ), at each vaccine frequency the 
steady state mean infectives were computed over ten net- 



work geometries. For the stochastic mean field (Fig 2(a) I, 
100 runs were done to extinction (Pj < 10~ 12 
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FIG. 1: Sample time series, a) Node fractions. Solid gray 
curve: infectives; dashed dark gray curve: susceptibles; solid 
black curve: vaccinated. '+' symbols indicates times of vac- 
cine pulses, b) Average degree by node class. Solid gray curve: 
infectives; dashed dark gray curve: susceptibles. p = 0.003, 
r = 0.002, q = 0.0002, A = 0.1, v = 0.0005, w = 0.04. 



means of the time series were computed. The mean field 
and full model are in excellent agreement. 

For the full system with rewiring (Fig. |2(d")j ), infectives 
were averaged over 8 x 10 4 MCS for a single adaptive net- 
work for the smaller vaccine frequencies (y < 0.00015), 
and error bars are the standard deviation. For larger 
v values the oscillations led to rapid die out, so the in- 
fectives were instead averaged over 100 runs computed 
to extinction, where the vaccine was turned on at time 
0. Error bars are the standard deviation of all the time 
points. This averaging includes transients but serves to 
illustrate the decreasing infective levels and large fluctua- 
tions due to oscillations as v is increased. The stochastic 

and 



mean field (Fig. 2(c) ) was computed as in Fig. 2 (a 



error bars are the standard deviations of the means. The 
vaccination frequency required to significantly lower the 
infected fraction in the adaptive network is much smaller 
than for the static network due to the interaction of vac- 
cination and rewiring. The mean field model accurately 
predicts the order of magnitude of vaccine required and 
the presence of oscillations. 

Figure [3] shows the dependence of the endemic state 
lifetime on the vaccine frequency for static (Erdos-Renyi) 
and adaptive networks. About two orders of magnitude 
less vaccine is needed in the adaptive case to significantly 
reduce the lifetime of the endemic state. For the full sys- 
tem, each point was obtained by averaging 100 runs for 
which the initial condition was the vaccine-free steady 
state and the vaccine was turned on at time zero. Mean 
field results were obtained similarly, but the time to ex- 
tinction was computed using a threshold for Pj of 10 -12 . 

To explain the efficacy of vaccination in the adaptive 
network, we examine the network structure in more de- 
tail, particularly the degree. Steady state degree distri- 
butions for an adaptive network with and without vac- 
cination are shown in Figure [D Results without vacci- 
nation were obtained similarly to those in 0. Results 
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FIG. 2: Average infectives versus vaccine frequency v for q — 0.0002, r — 0.002, p — 0.003, A — 0.1. |(a)| Mean field, w = 0. Gray 
and r ed c urves, respectively, are stable and unstable bra nch es of the deterministic mean field. Symbols are stochastic mean 
field. |(b) | Full system, w — 0. Curve is to guide the eyes, (c) Mean field, w — 0.04. Blue curve: stable steady state, red curve: 
unstable steady state, bla ck cu rve: stable periodic orbit, green curve: unstable periodic orbit for deterministic model. Symbols 
are stochastic mean field. 



(d) Full system, w = 0.04. (Color online. 
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FIG. 3: Dependence of endemic state average lifetime on vac- 
cine frequency v. Black open circles and dashed gray curve 
are full system and mean field, respectively, with rewiring 
(w = 0.04). Black closed circles and solid gray curve are 
full system and mean field, respectively, with no rewiring. 
p = 0.003, r = 0.002, q = 0.0002, A = 0.1. 



with vaccination were computed likewise but averaging 
over 9 runs to obtain better statistics for the vaccinated 
nodes which are present at very low levels (0.6% of the 
nodes). As shown in Fig. @Ia), susceptibles in the adap- 
tive network have a higher average degree than infectives 
due to rewiring. This is also apparent at the beginning 
of the degree time series for susceptibles and infectives in 
Fig. Dlb). Because the susceptibles typically have higher 
degree, random vaccination of susceptible nodes favors 
the higher degree nodes in the network and is therefore 
expected to be effective, as in previous studies of targeted 
vaccination [2^, [29|, [3(| • For a static network, in contrast, 



the high degree nodes are most likely to be infected and 
thus will rarely be selected for vaccination. 

The vaccination level used in Fig. @Jb) is small enough 
that it only slightly lowers the number of infected nodes. 
However, when vaccination occurs the rewiring becomes 
more effective because of the presence of "safe" nodes 
to wire to. The average infective degree is substantially 
decreased, and the susceptible degree is substantially in- 
creased. The vaccinated nodes have a broad degree dis- 
tribution with a very high average degree. For the pa- 
rameter values in Figure HJb) , vaccinated nodes partic- 
ipate in links as frequently as susceptibles, even though 
there are two orders of magnitude fewer vaccinated than 
susceptibles. 

In this model, we randomly vaccinated a fixed frac- 
tion A of the susceptible nodes during each vaccination 
event. Information about the network structure was not 
required. This corresponds to the case where vaccine is 
inexpensive and anyone who is not infected can be vac- 
cinated. If the vaccine were in limited supply, it would 
be desirable to use some knowledge about the social net- 
work structure (such as the degrees) to target the vaccine 
more carefully. In that case, the vaccine pulses might be 
a fixed number of susceptibles rather than a fixed frac- 
tion. However, further study is needed to determine the 
optimal vaccination strategy. Another extension would 
be to relax the assumption that nodes have full knowl- 
edge of others' disease status 0, H|. 

In summary, we studied the effect of Poisson vaccina- 
tion on epidemic spread in an adaptive network. Vacci- 
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nation was far more effective (by two orders of magni- 
tude for the parameters chosen) in an adaptive network 
than a static one due to the interaction of vaccination 
and rewiring. The network adaptation led to a higher 
average degree for susceptible nodes, and random vacci- 
nation of susceptibles exploited this heterogeneity. The 
adaptive rewiring then became more effective because of 
the presence of "safe" vaccinated nodes to rewire to. 
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FIG. 4: Degree distributions from Monte Carlo simulation for 
p = 0.002, r = 0.002, w = 0.04. (a) No vaccination, (b) With 
vaccination: v = 0.00015, A = 0.1, q = 0.0002. Degree distri- 
bution of vaccinated nodes is very broad, extending beyond 
the figure domain, with an average degree of 579. Solid gray: 
infectives; dashed: susceptibles; solid black: vaccinated. 
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